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ABSTRACT 

We present new UBVI C CCD photometry of the young open clusters Trum- 
pler 14 (Tr 14) and Trumpler 16 (Tr 16) in the r] Carina nebula. We also identify 
the optical counterpart of Chandra X-ray sources and Two Micron All Sky Survey 
point sources. The members of the clusters were selected from the proper mo- 
tion study, spectral types, reddening characteristics, and X-ray or near-IR excess 
emission. An abnormal reddening law Ry, c i = 4.4 ±0.2 was obtained for the stars 
in the r\ Carina nebula using the 141 early-type stars with high proper motion 
membership probability (P M ^ 70%). We determined the distance to each cluster 
and conclude that Tr 14 and Tr 16 have practically the same distance modulus 
of V Q - M v = 12.3 ± 0.2 mag (d = 2.9 ± 0.3 kpc). The slope of the initial mass 
function was determined to be T = -1.3 ± 0.1 for Tr 14, T = -1.3 ± 0.1 for Tr 
16, and T = —1.4 ± 0.1 for all members in the observed region for the stars with 
logm ^ 0.2. We also estimated the age of the clusters to be about 1-3 Myr from 
the evolutionary stage of evolved stars and low-mass pre-main-sequence stars. 

Subject headings: open clusters and associations: individual (Trumpler 14, Trumpler 
16) — stars:luminosity function, mass function — stars: pre-main sequence 
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Introduction 



The universalit y of the initial mass function (IMF) has been debated from the time 



that ISalpeterl (ll955f ) first introduced the idea of an IMF. The massive part of the IMF can 
be represented by its slope from the inflection point at ~ 1M Q to the most massive stars and 
is described by its slope (T = ^/°g ^ ) , which is the so-called Salpeter-type IMF. Several star 
forming regions show somewhat shallower slopes of T in the massive part, particularly for 
the most massive star forming clusters, even though we take the observational uncertainties 
in the photomet ric calibrations into a ccount. For example, the Arches cluster shows 



r — —1.1 ± 0.2 ( Espinoza et al 



( ISung &: Bessell 



20091 ) and the core of NGC 3603 shows T = -0.9 ± 0.1 



2004 ). However, the main reason for the difficulty in investigating the IMF 



of these clusters is that they are either highly reddened or very distant. 

The 7] Carina nebula (NGC 3372) is the largest nebula in the southern sky and is one 
of the most interesting regions in the Galaxy. More than 60 O-type stars and several young 
open clusters are located in or near the nebula, including Tr 14 and Tr 16 in the bright part 
of the nebula. Moreover, rj Carinae, which is known as one of the most massive stars in the 
Galaxy, is near the center of Tr 16. Although the Arches cluster and NGC 3603 appear to 
be more massive, the young open clusters in the rj Carina nebula offer many advantages. 
They are relatively close and less reddened, and are therefore able t o be inve s tigated in 



more det a il. M a ny previous investigators ha v e studied these clust ers 



Walborn 



(120021 ) . iLevato k Malarodal (119821 ) 



Morrell et a 



( 119931 ) presented spectral types of the stars in the region. 



Walborn ( 



1988), and 



1973 



1995 



Massev fc Johnson 



Cudworth et al. 



( 119931 ) performed 



mag. Previous studies i nclude optical 



1993 



1987 



Tapia et al. 


2003 


Tapia et al. 


2003 



Carraro et al. 



Ascenso et al. 



2007 



probability of s 


tars down to V = 16 


(1 Vazquez et al. 


1996: 


Cudwon 


;h et al. 



2004) , near-IR (NIR) JHK photometry (ISmith 



Preibisch et al. 



20111 ). or X-ray observations 
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( lAlbacete-Colombo et al.l 120081 ). iMassey fc Johnson! ( 
(r = -1.3 ± 0.2) for Tr 14 and Tr 16 for M ^ 6M . 



993) obtained t 



r e slop e of the IMF 



Smith &: Brooks 



( 120081 ) summarized 



the obser vational evidences of o ngoing; star formation activity in the rj Carina nebula. 



Recently, 



Townsley et al. 



(120111 ) (the Chandra Carina Complex Project - CCCP) published 
a vast amount of information on the stellar content and star formation activity in the rj 
Carina nebula region. 

As these clusters are in the line of sight toward the tangential point of the Sagittarius 
spiral arm, it has been questioned as to whether the young open clusters in the 77 Carina 
nebula are at the same distance or not. Tr 14 and Tr 16, in particular, have been at 



the center of this controversial question. In addition, it is also a r 



clusters ex. 



libit an abnorma 



The 



et al. 



1996 



1980 



Smith 



Tapia et al. 



2003 



1987 



reddening law. (Feinstein et al. 



Tap ia et 



Carraro et al] 



al. 



1988 



2004 : 



1973 



Herbst 



Massev &: Johnson 



Ascenso et al. 



gued that t 



1993 



rese two 



1976 



Forte 



1978 



Vazquez et al 



20071 ). Most pho tometri c stud ies 



Smith! ()2006j) 



have estimated the distances of these clusters to be 2.0-4.0 kpc. Recently, 
obtained 2.35 ± 0.5 kpc for the distance to 77 Carinae using the proper motion of the 
Homunculus Nebula around 77 Carinae. 

The aims of this study are (1) to judge whether the reddening law in the 77 Carina 
nebula is abnormal or not, (2) to determine the distance to the young open clusters in 
the nebula and decide whether they are at the same distance or not, and (3) to obtain 
the IMF of the young open clusters in the nebula down to the inflection point (~ 1M Q ) 
and determine its slope. In Section 2 we present our new optical data and describe the 
data sets used in this work. In Section 3, we determine the reddening law and distance 
modulus. The membership selection criteria are presented in Section 4 where we construct 
the Herzsprung-Russell (H-R) diagram and determine the age and masses of individual 
stars in the clusters. The IMF is derived in the same section. The summary is presented in 
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Section 5. 



2. Observation 
2.1. Optical data 

UBVI C observations were performed on 1997 March 2, June 23, 1999 February 5 and 
6 using the lm telescope at Siding Spring Observatory with a SITe 2K CCD camera. The 
field of view was 20 '.5 x 20 '.5. We obtained two short and four long exposure images for 
each filter. The average seeing was ~ 2" during all observing runs. The total observed area 
covered 22'A x 20 f 8 of the nebula. 



We performed point-spread-function (PSF) photometry using the DAOPHOT package 
in IRAF o An aperture correction was applied to produce the equivalent magnitude for a 
7" radius. We tran sformed our observed data to the SAAO standard system using SAAO 



E5 and E7 regions (ISung fc Bessell 
presented in Table [U 



2000). The photometric data down to V = 19 mag are 



EDITOR: PLACE TABLE □ HERE. 



We compared our p hotometry wit 



Table [2l Three studies (IVazquez et al 



i previous studies, and t he results are presen ted in 



based on CCD PSF photometry while 



1996 



Tapia et al. 



Massey fc Johnson 



2003 



Carraro et al. 



20041 ) were 



(119931 ) performed simple aperture 



photometry. Others were based on photoelectric photometry. In the comparison, we have 



Tmage Reduction and Analysis Facility is developed and distributed by the National 
Optical Astronomy Observatories which is operated by the Association of Universities for 
Research in Astronomy under operative agreement with the National Science Foundation. 
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excluded stars that deviated by more than 2.5cr from the mean to avoid the inclusion 
of variables or optical doubles. Our photometry is in good agreement w ith previous 



Vazquez et al. 



fll996[ ). But although 



photoelectric photometry and CCD photometry by 
we find no significant difference in the photometric zero point in F, we find a large scatter 
in the comparison with t he photoelectric data and simple aperture photometric data of 



1 



Massey fc Johnson! (119931 ). Such a large scatter may be caused by difficulty in correcting for 



the sky background due to spatially varying nebulosity. There are no s ystematic difference s 



in c olor between our pho tometry and previous photometry except with 



and 



Carraro et al 



fl2004h . 



Tapia et al 



fl2003f ) 



The comparison with previous CCD photometry is shown in Fig ure [U The top pan el 



shows the comparison with the simple aperture photometric data of 



Massey fc Johnson 



(119931 ). This shows the largest scatter, as mentioned aboye. T 



comparison with the CCD photometry of 



Vazquez et al 



re second panel shows the 



(119961 ). There is no significant 



systematic difference with their data except in (V — I) for some stars around Tr 14. Most 
of these stars are brighter than V < 14.5 mag and located in a small area (—3^5 < A R.A. 



< -7'.5, 7'.0 < A decl. < 
with CCD photometry by 



0:5) . In the 



Tapia et al. 



;hird and the last panels we sh ow the comparison 



(120031 ) and 



Carraro et al. 



(120041 ) . Their data show 



large differences in the photometric zero points as well as a large scatter. 



EDITOR: PLACE TABLE El HERE. 



EDITOR: PLACE FIGURE Q] HERE. 
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2.2. X-Ray data and 2MASS data 



CCD coordinates have been transformed to the equator ial coordinate system using the 



2006|). A total 



Two Micron All Sky Survey (2MASS) point-source catalog (ISkrutskie et al.l 
of 4615 out of 6235 optical sources have been identified with 2MASS point sources within a 
search radius of 1". Figure [2] shows the (J — H) versus (H — Ks) color-color diagram for 



the stars with e( J — H : H — Ks) (= J tj 2 + 2e# 2 + e^ s 2 ) < 0.1 mag. The ma in-sequence 



(MS) relation of ISung et al 



( 120081 ) and the reddening law of iFitzpatrickl ( 119991 ) have been 



adopted and used in the figure. Thirty-four stars were found to have excess emission in 
(H — Ks) color and are classified as NIR excess stars. 



EDITOR: PLACE FIGURE M HERE. 



X-ray observations are the most powerful tool for selecting P MS members (jSung et al 



2004). The Chandra X-ray point source catalog publish ed by 



(120081 ) and the CCCP X-ray source catalog published by 



selecting the X-ray emission stars. 



Albacete-Co 



Broos et al 



Albacete-Colombo et al. 



ombo et al. 



(1201 ll ) were used in 



(120081 ) also published X-ray 
sources in the 17' x 17' area centered on R.A.= 10 h 44 m 47593 and decl.= -59°43 m 54 f !21, 



which covers the entire region of Tr 16. The CCCP provi ded a catalog of 14369 X 



Albacete-Colombo et al 



ray 



(}2008|). 



sources in a 1.42 deg 2 area including the region observed by 
There were 3338 CCCP sources in our field of view. Both X-ray catalogs supplied the 
2MASS counterpart of X-ray sources. First, we used th eir identification and found 964 
and 373 X-ray emission stars in the CCCP catalog and I Albacete-Colombo et al.l (120081 ) 



respectively. Next, we searched for optical counterparts of X-ray sources without a 2MASS 
ID using a search radius of l'/2. If two or more stars were found within the radius, we 
assigned the closest star as the optical counterpart of the X-ray s ource. We identified 192 



and 96 stars as X-ray emission stars from the CCCP catalog and lAlbacete-Colombo et al. 
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( 120081 ). respectively. As a result, a total of 1184 stars were selected as X-ray emission stars. 



Among them, 441 stars were identified as X-ray sources from both X-ray catalogs. 



3. Reddening and Distance 



3.1. Photometric Diagrams 

Figure |3] shows the color-magnitude diagrams (CMDs) of stars with e < 0.1 mag in 
the observed region. Thi c k line s represent the reddened zero-age main-sequence (ZAMS) 



relation from 



Sung et al 



(119971 ) for E(B - V) = 0.61 mag. The stars brighter than V = 9 



mag were saturated in our observations. We therefore used the mean magnitudes and colors 
from the SIMBADy database for these stars and they are marked as open circles in the 
figures. 



T he stars with hig 



study ( iCudworth et al. 



i me mbership probability (P M ^ 70%) from the proper motion 



19931 ) were classified as proper motion members (see § I4.ip and are 



marked as filled circles. X-ray emission stars and NIR excess stars are marked as crosses 
and triangles, respectively. Most of the proper motion members are scattered around the 
reddened ZAMS line. The relatively large scatter implies a wide range of reddening among 
the cluster members. There are several red stars (V — I > 1.0 and V < 3.75(V — I) + 8.0) 
with a high proper motion membership probability. Seven stars (ID 1817, 2599, 2639, 2754, 
2830, 3337, and 3911) seem to be foreground late-type stars with similar proper motion 
vectors by chance to the clusters. Four stars (ID 2022, 2033, 4211, and 4238) are near 
the boundary of bright nebula; they are probably highly obscured (E(B — V) > 0.9 mag) 
early-type stars. Six additional stars (ID 677, 849, 1027, 1374, 2462, and 4044) have similar 
photometric characteristics as the four stars mentioned above. The star ID 849 is an optical 



2 http: / / simbad.u-strassbg.fr/simbad/ 



-9 - 



double, and shows a variability of AV ~ 0.25 mag from four observations. In addition, the 
star ID 2572 (V = 16.265, V - I = 1.375, B - V = 0.461, U - B = -0.440) shows very 
strong UV excess with a NIR excess; this star appears to be a classical T Tauri star with 
active accretion. 

Most X-ray emission stars are fainter than V = 15 mag and redder than the reddened 
ZAMS line by about 0.6 mag. They constitute a well-defined sequence of PMS stars in 
the clusters. However, a non-negligible number of X-ray emission stars can also be found 
around the reddened ZAMS line in the (V, V — I) diagram. They seem to be X-ray active 
field stars in front of the rj Carina nebula because their V — I color is not affected by 
the UV excess and late-type stars retain their X-ray activity for a long time. There are 
several X-ray emission stars redder than most PMS stars with X-ray emission; two of 
them (ID 2022 and 4211), as already mentioned, are highly reddened early- type stars with 
high proper motion membership probability. The star ID 2412 seems to be a foreground 
late-type star because its position in the (U — B,B — V) diagram is close to the unreddened 
ZAMS line. The X-ray emission star ID 1836 seems to be a late-type member with a strong 
UV excess. The other three stars (ID 3729, 4081, and 4212) are very difficult to justify as 
members from optical data alone; they seem to be highly reddened early-type stars in the 
(J — H, H — K s ) diagram and, in addition, are in the dark lane to the south of Tr 16. The 
membership of these presumably highly reddened stars, including the four stars mentioned 
above, is uncertain, and they are for the time being, not included in the member list. 

EDITOR: PLACE FIGURE [3] HERE. 

Figure H] shows the (U — B , B — V) color-color diagram. Early-type members without 
proper motion membership probability were selected from the locus of early-type members 
in Figure [5] (see Section 13.21) and are marked as open squares in Figure H] and [51 The 
reddening of proper motion members and bright members is between E(B — V) — 0.36 
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and 0.9. There are several stars with E(B — V) < EiB — V)f g = 0.36 mag and with high 
proper motion membership probability (P^ ^ 70%) including six early-type stars (ID 2154, 
2377, 2475, 2693, 3223, and 4156). Were they actual members of the r] Carina nebula, 
they should be on the blue part of the reddened ZAMS whereas they lie in the middle of 
the MS band in the (V, U — B) CMD. This implies that they are more likely foreground 
stars in the spiral arm between t] Carinae and the Sun (see § 13. 2p . There are also several 
proper motion members in the red part of the E(B — V) = 0.9 line. Blue stars with 
Q = (U — B) — 0.72(5 — V) ^ —0.5 mag are highly reddened early-type stars, while most 
red stars near the E(B — V) — 0.9 line are foreground late-type stars. In addition, there are 
many X-ray emission stars between ( B — V) = 0.9 - 1.3 mag and (U — B) = 0.0 - 0.8 mag. 



These are PMS stars with UV excess (ISung et al. 
16 are actively accreting. 



19971 ) implying that the PMS stars in Tr 



EDITOR: PLACE FIGURE 1 HERE. 



3.2. Reddening law 



It has long been debated whether or not Tr 14 and Tr 16 are at the same distance. The 
reddening law may play a key role in this debate. The total extinction in the V band can be 
determined using the relation Ay = Ry x E(B — V). The so-called normal reddening law of 
Rv = 3.1 ± 0.2 is cons istently obtained for stars in the solar neighborhood ( jGuetter fc Vrba 



1989 



Lim et al 



20111 ). while many previous studies noted that the total-to-selective 
extinction ratio -Ry toward the n Carina region is anomalously higher 



Hcrbst 



Smith 



1976 



20021 ) 



brte 



1978 



The et al 



Vazquez et al. 



1980 



Smith 



1987 



Tapia et al 



1988 



(Feinstein et al. 


1973 


; Vazauez et al. 


1996 



( 119961 ) used the membership probability (P^) flCudworth et al. 
19931 ) as the membership selection criterion and determined E(B — V)f g = 0.33 mag and 
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R 



V,cl 



also 



4.7 ±04 for Tr 14 [A v = 3.1 x E(B - V) fg + R v ,d x [£(5 — V) — E(B - V) Jg \, see 



Forte! fll978h ]. 



We have determined the reddening law toward the 77 Carina nebula using t he relation 



between Rv and color excess ratios, R\ 



2.45 eZ-v) flGuetter & Vrba 



19891 ) . The 



£"(5 — V) reddening of individual early-type stars is determined in the (U — B , B — V) 
diagram using the slope of the reddening vector fpjz^ = 0.72. The reddening in (V — I), 
E(V — I), is the difference between the observed and the intrinsic (V — I). The intrinsic 



(V — I) is d etermined from t 



Appendix of 



re re lation between (B — V) and (V — I)o for MS stars (see 



Sung &: Besselll (Il999h ) . We selected early-type stars to be those with V < 15.5 



mag and Q = {U - B) - 0.72(fi - V) S -0.5 mag or with V < 15.5 mag, (U - B) < -0.1 
mag and (B — V) < 0.45 mag. The early-type membership selection criteria are P M ^ 70% 
and 0.36 ^ E(B — V) ^ 0.9 mag. Stars with E(B — V) > 0.9 mag were excluded to avoid 
contamination by PMS stars with strong UV excess or background early-type stars. A 
total of 141 early-type members were selected and used to determine the reddening law and 
distance modulus of the r] Carina nebula. Figure |5] shows the relation between E(B — V) 
and E(V — I). We determined the color excess ratio due to the intracluster dust to be 

1.80 ± 0.10 and E(B - V) fg = 0.36 ± 0.04 mag using le ast 



E(V-IU 
E{B-V) cl 



E(V-I)-E(V-I) fg ■ 
E{B-V)-E(B-V) fg - 



squar es. The color excess ratio e{b-v) ^ or a norma l reddening law is 1.25 ( Dean et al. 



19781 ). The total-to-selective extinction ratio in the 7] Carina region (Rv,d) is, therefore, 
4.4 ±0.2. 

If the six early-type stars with E(B — V) < 0.36 mag are included, the resultant 
foreground reddening E(B — V)f g and reddening law Ry are 0.25 ± 0.09 and 4.0 ± 0.25, 
respectively. Although we could derive the same distance modulus in the case of a smaller 
Ry, the reddened ZAMS lines are not well fitted to the distribution of early- type members 
in Figure [3j This fact supports the suggestion that these less reddened early-type stars 
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are foreground stars in the Sagittarius spiral arm between Vq — My = 11.3 -11.8 mag (see 
Section 13TTj) . 



Vazquez et al. 



(119961 ) also applied the same method to determine the reddening law of 
the nebula. Although there are small differences in the membership selection criteria, their 
result agrees well with ours. The small difference between our result and theirs appears to 



be caused b y the number o 



were used in 



Vazquez et al. 



stars used in the determination of the reddening law (53 stars 
(jl996[ ). while 141 stars were used in this work). 



EDITOR: PLACE FIGURE HERE. 



As there is no reliable photometric way to estimate the reddening of PMS stars, we 
have to rely on the reddening map derived using the early-type stars. Figure [6] shows 
the reddening map. The dashed line represents E(B — V) — 0.5 mag and the solid lines 
represent E(B — V) = 0.6 and 0.7 mag, respectively, from thin to thick lin e. The spatial 



variat ion of E(B — V) in Figure E] is well matched by th e reddening map of 



(119731 ) and the extinction map derived by 



Feinstein et al 



Smith! (119871 ). The highly reddened region is 



to the west of Tr 14 (AR.A.> -lO'.O and 4f0 < Adecl.< 8(0, where AR.A. and Adecl. 
represent the angular di stance in R.A . and D ec, from rj Carinae) and is well matched by 



the CO emission map of 



Brooks et al 



(I2003h . 



EDITOR: PLACE FIGURE E HERE. 



3.3. Distance 



In order to derive the distance modulus of each cluster, we should first determine the 
membership containing radius of each cluster. However, as it is still uncertain whether Tr 



13 



14, Tr 16, and Cr 232 are dynamically independent clusters, it is therefore very difficult 
to set the boundaries of each cluster. We decided to set the centers and boundaries of 
each cluster so as to include most of proper motion members (P M ^ 70%). The center and 
radius of each cluster are A R.A.= — 8fl, A decl.= 7'2, and radius = 4'.2 for Tr 14 and A 
R.A.= -i:5, A decl.= -2'.0, and radius = 6'.0 for Tr 16 (see Figure 0). 



EDITOR: PLACE FIGURE [7] HERE. 



Although the proper motion study and spectral classification studies can supply 
good membership criteria, it is still debate d whether Tr 14 and Tr 16 are at the same 



distance or not. Two rec ent optical stu d ies ( 



at different conclusions. 



Carraro et al. 



apia et al. 



2003 



Carraro et al. 



2004 ) arrived 



(120041 ) determined the reddening law in three 
different ways and adopted an average value. They adopted Ry = 3.48 and 4.16 for Tr 
14 and Tr 16, respectively, and obtained a different distance modulus for each cluster 
(Vo — My = 12.3 ± 0.2 mag for Tr 14 and 13.0 ± 0.3 mag for Tr 16). However, they did 
not take into account the difference, in star-forming regions, of the reddening law between 
the general interstellar medium and the intracluster medium. In addition, the number of 
members used in the determination of the re ddening law was v ery small (10 stars for Tr 14 



Tapia et al. 



and 14 stars for Tr 16). On the other hand, 
mag using Ay = 1.39E(V — J) to avoid the anomalous reddening 
and concluded that Tr 14 and Tr 16 are at the same distance. 



(120031 ) obtained V - My = 12.1 
aw at optical wavelen gths 



Vazquez et al. 



(1996) 



determined Vq — My = 12.5 ± 0.2 mag for T r 14 by apply ing the abnormal reddening law 



they determined (see Section |3~2]) . Recently, ISmithl ( 120061 ) determined the distance to 77 



Carinae as d — 2.35 ± 0.05 kpc (Vo — My = 11.85 ± 0.05 mag) from the proper motion of 
the Homunculus Nebula ejected from 77 Carinae. 



EDITOR: PLACE FIGURE M HERE. 
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We calculated the total extinction in the V band of individual early-type memb ers 



19971 ). and 



using the reddening-corrected colors and the adopted ZAMS relation ( jSung et al. 
constructed the distribution of distance moduli for each cluster. We find the same distance 
modulus for Tr 14 and Tr 16 (Vo — My = 12.3±0.2 mag, i.e. d = 2.9±0.3 kpc) and conclude 
that Tr 14 and Tr 16 are at the same distance within the observational errors. Figure |8] 
shows the reddening corrected CMDs superimposed with the ZAMS relation shifted by the 
adopted distance modulus. Although it is very difficult to estimate the error in the ZAMS 
fitting, we can expect an error of at least 0.1 mag. In addition, the ZAMS relation itself 
may be uncertain by about 0.1 mag. The error in the distance modulus determination 
is, therefore , abou t 0.2 mag. Our dist a nce m odulus is in good agreement with that by 



apia et al. 



( 120031 ) and IVazquez et all ( 119961 ) , but disagrees with that by 



Carraro et al. 



( 120041 ) who estimated different distances to Tr 14 and Tr 16. This difference is caused by 



the applied reddening correction, i.e. the difference in Ay of individual early-type stars. 

There is a non-negligible difference in the derived distance between the two methods, 
the proper motion of the expanding nebula ejected fro m 77 Carinae (the Homu n culus Neb ula) 



and th e ph otomet r ic analysis of stars in the nebula. 

~~ith| 



(119991 ) . and 



Smith! (12001 



Allen fc Hillierl (119931 ) 



Meaburn 



20061 ) performed spectroscopic analyses of the Homunculus Nebula 
and determined 2.2 ± 0.2 kpc, 2.3 ± 0.3 kpc, and 2.35 ± 0.05 kpc, respectively. These are 
quite consistent each other within th e expe cted error. On the othe r hand, the distance 



moduli determined by 



Vazquez et al. 



Tapia et al 



(120031 ) are comparable with 



(11999 ) and 

our result of d = 2.9 ± 0.3 kpc. Because these latter two studies and our work incorporate 
the abnormal reddening law toward the 77 Carina nebula, the consistency of the distance 
moduli from the photometric studies suggests that our distance determination is realistic. 
Furthermore, the possibility cannot be ruled out that the distance from the proper motion 
studies of the ejected nebula might be affected by the assumption of a constant expanding 
velocity and symmetrical shape of the Homunculus Nebula. If the expansion velocity of the 
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ejected material varied after the initial ejection, or if the Homunculus Nebula is asymmetric 
in shape, the proper motion of the expanding nebula may be under estimated. In addition, 



1997|), the influence of the 



as rj Carinae is known to be a binary system ( IDamineli et al. 
companion star on the expansion of the Homunculus Nebula should be taken into account. 
Therefore, although our distance modulus is somewhat larger than the distance determined 
from the proper motion of the ejected nebula, our value reasonably explains the location of 
the MS band. 



Recently 



Artigau et al 



(1201 ll ) studied the proper motion of the Homunculus Nebula 
and determined the epoch of the eruption of rj Carinae. They assumed the distance to t] 
Carinae as 2.3 kpc. Although most of parameters (projected velocities, orientation angles, 
and projected distances from rj Carinae) may be changed if we adopt the distance to r\ 
Carinae to be 2.9 kpc, the epoch of eruption is nearly the same because the projected 



velocity of the nebula generally dominates the radial velocity. 



4. H-R diagram 

4.1. Membership selection 

The reddening law was used as the primary membership selection criterion for 
early-type stars. In Figure most of the early-type stars with P M ^ 70% lie along 
the line of Rv,ci — 4.4. The lower boundary of membership selection was set to be 
AE(B — V) — 0.06 mag at a given E(V — I). We classified stars above the lower boundary 
and with 0.36 ^ E(B — V) ^ 0.9 mag as early-type members. Early-type members 
selected using the reddening law and without proper motion membership probability are 
marked as open squares in Figure [H From this criterion, some stars with low membership 
probability (20% ^ P M < 70%) were also selected as members because nearly face-on binary 
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systems may have a low proper motion membership probability. For instance, HD 93250 
(03.5V((f+)), = 24%), one of the earliest type stars in the y Carina region, was selected 



as a m ember. The star shows unusually strong X-ray activity (jEvans et al 



2004 : 



Naze et al. 



20 111 ) and is a suspected binary syste m. But it is still n ot confirmed as a binary star due 



to a lack of radial velocity variations (iRauw et al. 
using the reddening criterion, the total extinction Ay was ca 
A v = 3.1 x E(B - V) fg + R v ,d x [E(B -V)- E(B - V) fg ) ((Forte 



20091 ) . For early-type members selected 
culated ind ividually using 



19781 ) 



For intermediate or low- mass PMS stars, E(B — V) was estimated from the reddening 
map of Figure O Stars with reddening-corrected colors within |A(V — J)o| ^ 0.15 mag 
and \A(B — V)o\ S 0.15 mag from the ZAMS relation, and with membership selection 
criteria (X-ray emission, P M ^ 70% or NIR excess) were classified as members and stars 
with 20% ^ P M < 70% were classified as candidates. To select the PMS members, we se t 



the PMS locus in Figure [3] to include most of the X-ray emission stars (jSung et al. 



20081 ) 



The thick dashed line in Figure [3] represents the PMS locus. We classified stars satisfying 
the above membership selection criteria (i.e. X-ray emission, P M ^ 70% or NIR excess and 
falling within the PMS locus) as members. The stars in the PMS locus without membership 
criteria were classified as PMS candidates. The spatial distribution of members is shown in 
Figure [3 



4.2. Color— T e fj Relation and Bolometric Correction 



In order to construct the H-R diagram, we transformed the reddening-corrected 
colors, spectral type, and My to effective temperature and bolometric magnit ude. The 



tempe rature of the 02 supergiant HD 93129A (02If*) is somewhat uncer tain. 



1996) determined the temperature of HD 93129A to be 50,500 K, while 



Puis et a. 



Repolust et al. 



(120041 ) estimated it to be 42,500 K. However, a temperature of 42,500 K is much too low 
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for spectral type 02. As there are no confirmed 02If* stars except HP 93129A in the 
Galax y, we referred to the temperature of 02If* stars in the LMC 



Doran fe Crowther 



( 1201 ll ) determined the effective temperature of two 02If* to be about 50,000 K using N IV 
and N V lines, instead of He lines. We, therefore, assumed T e ff = 50,000 K for HD 93129A. 
For 03 to 08 s tars, we applied the s pectral- type-temperature relation and the bolometric 



correction from iMartins et al. 



(120051 ) . For 09 to (V — J)o < 1.6 mag stars, we adopted 



the co lor-temperature and temperature-bolometric correction relations from 



Bessell et al. 



( 119981 ). For (V — I)o < mag stars, (B — V)o and {U — B) were used as the temperature 
indicators. For ^ (V — I) 5= 1.6 mag stars, (V — I) was used as the temperature 
indicator. For (V — I)q > 1.6 mag stars, w e have a dopte d the empirical color-temperatur e 



relation and bolometric correction scale of 



Besselll (119911 ). bavidson k Humphrevsl (119971 1 



estimated the total luminosity of 77 Carinae (LBV) to be ~ 1O 6 ' 7 L if the star is at 2.3 
kpc. For 7] Carinae, we adopted the effective temperature of log (T e yj) = 4.128 from the 
color-temperature relation and the luminosity of 8 x 1O 6 L for the case of d = 2.9 kpc. 
For HD 93162 (WN7), we _ adopted the effective temperature of log (T e ff) = 4.491 and the 

( 11995k 



luminosity of 10 L & from 



Crowther et al 



EDITOR: PLACE FIGURE M HERE. 



4.3. The Initial Mass Function 



The H-R diagram with several evolutionary tracks is shown in Figure [9l In the 
figure, most stars with log(T e ^j) > 4.0 follow a well-defined sequence near the ZAMS. To 



determine the mass of individual sta rs, we adoptee 



the non-rotating evolution models of 



Schaller et al. 



the initial rotational velocity of 330 km s 1 of 



three different evolution models 



(1992J), ( n) the rotating models with 



Brott et al. 



(1201 lh for 5-60M© and the 
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of 



Ekstrom et al. 



assumed 75 M ( IGamen et al. 



no n-rotating mod e ls of ISchaller et all ( 119921 ) for the other stars, and (iii) the rotating models 



(20 l l|). Because th e mass of LBV and WN7 stars are very uncertain, we 



2008|) and 100 M Q for HD 93162 (WN7) and rj Carinae, 
respectively. Although both stars seem to be binary systems, we only considered the masses 
of the primaries because photometry alone cannot take into account the secondary of a 
binary system. We also assumed the upper mass limit to be 100 M & . There are two more 
stars [HD 93129A (02If) and HD 93250 (03.5V((f+)))] with even higher mass. These two 
stars are brighter than the evolutionary track of a 120 M & star and suggests that the mass 
of the evolved stars (H D 93162 and y Carinae ) may be more massive than the upper mass 



limit (100 M Q , see also 



Crowther et al. 



fcoioh ). 



EDITOR: PLACE FIGURE [TO] HERE. 



We calculated the IMF £(= N/A logm/area) for the case of A logm = 0.2. The error 
bars are based on Poisson noise. To minimize the binning effect, we also calculated the 
IMF by shifting 0.1 in logm and then calculated the slope of the IMF T(= d log £/d log m) . 
The slope of the IMF (T) calculated for members and candidates —1.3 ± 0.1 (standard 
error) for Tr 14, —1.3 ± 0.1 for Tr 16, and —1.4 ± 0.1 for all observed regions down to our 
observational limit (logm ^ 0.2, Table |3]case (1)). 

Figure PTOl shows the IMF of Tr 14, Tr 16, and all observed area. In the figure, the 
IMF calculated for members and candidates is marked as filled circles and that calculated 
for members-only is marked as open circles. All IMFs are the same as for Tr 14. As Tr 14 
is a compact cluster and most members are assembled in a small area, we can expect the 
contamination of field stars to be very low. On the other hand, the IMF of members and 
candidates of Tr 16 is slightly steeper than the IMF calculated for the members-only. The 
difference is only apparent for logm < 0.6. From the figure, we can see that the IMF for 
members and candidates and the IMF for members-only differ somewhat for logm < 0.6. 
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That is caused by the observational limit of the proper motion study (ICudworth et al. 



19931 ) . Therefore, the IMF calculated for members including the candidates should be 



considered as an upper limit of the IMF of the rj Carina nebula region (Table [3] case (1)) 
and the IMF of the members-only can be considered as a lower limit (Table [3] case (2)). 

If all stars in the PMS locus, regardless of membership selection criteria, are used in 
the calculation, the IMF may be overestimated due to field star contamination. To subtract 
the contribution of field stars, we should choose a control region representative of the field 
population toward the rj Carina nebula. However, the choice of the field region is very 
difficult, especially for open clusters in the galactic plane. If we choose a region far from 
the nebula, the background population will be counted more than the nebula region. On 
the other hand, the region near the i] Carina nebula may represent both foreground and 
background populations, but faint halo members of the cluster may also be subtracted as 
part of the "foreground population" . 

In order to subtract the field contribution and compare the IMF before and after the 
subtraction, we chose a small region in the N W part of the observed region. The region is 
shown in Figure [7] (A R.A.< — 13'5, A decl.> 0-5). The l ocation of the reg i on co incides 



with the giant molecular cloud Car I ( iBrooks et al. 



20031 ). ISmith k Brooksl (120081 ) showed 



that the bright arc in Ha is well matched with the curved edges of both 13 CO and 3.3 /im 
PAH emission. These facts allow us to expect that star formation has not yet propagated 
into this region, and we can expect that no members of the rj Carina nebula exist in the 
region. In addition, as the giant molecular cloud is very close to the rj Carina nebula, stars 
found in this region should well represent the population of foreground stars. 

We re-calculated the slope of the IMF, T, after subtraction of the field contribution 
down to our observational limit (Table |3] case (3)). However, as there was no significant 
difference from the slopes of the IMF in Table [3] we conclude that the slope of IMF, 
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T, is practically the same within the errors, with and without subtraction of the field 
contribution. The field contamination in the observed region has no serious influence on the 
slope of the IMF. 



The IMF of the rj Carina nebula is very similar to that of NGC 6231 (r = — 1.2 ± 0.4; 



Sung et al. 



19981 ). but is slightly shallower than t he IMF of the sol ar neighborhood, the 



Pleiades, or NGC 2264 (r « -1.7; see Figure 6 of 



Sung et al. 



fl2004|)). 



1: 



4.4. Age and Age Spread 



The presence of the LBV (77 Carinae), WN7 (HD 93162), and several 03 stars indicates 
that both Tr 14 and Tr 16 are very young clusters . Most previou s estim ates of age and age 



spread were based on stellar evolutio nary models. 



for the most massive stars in Tr 14. 



Vazquez et al. 



Degioia-Eastwood et al. 



9961 ) estimated 1.5 Myr 



( 1200 ll ) estimated the age of 



both Tr 14 and Tr 16 to be 2 - 3 M yr with an age spre ad of 5 - 6 Myr for massive stars 



and 10 Myr for low-mass PMS stars. iTapia et al. 



Carraro et al. 



jgOOa ) estimated 5 Myr for Tr 14 and 3 



(12004f l estimated 2 Myr for Tr 14 and 5 Myr 



6 Myr for Tr 16. Most recently, 
for Tr 16. 

In Figure El tw o isochrones interpolated from the stellar evolution models of 



Schaller et al. 



dl992| ) are superimposed. HD 93128 (03.5V((f+))) and HD 93129A (02If*) 
in Tr 14 and HD 93205 (03.5V((f+))), HD 303308 (04V((f+))) and ID 2022 (=Y398, 
03-04If) in Tr 16 are very close to the ZAMS and fit well to the 1 Myr isochrone, while HD 



93162(3 (WN7) in Tr 16 is well fitted to the 2.5 Myr isochrone. Therefore, the age of massive 



For HD 93162, 



Hamann et al 



(120061 ) obtained an effective temperature of 50,000K 



an d luminosity o f 10 6 55 L^ for the case of Vq — My = 12.55. If the effective temperature 



of 



Hamann et al. 



( 120061 ) and the luminosity corrected for the case of Vq — My = 12.3 are 
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MS stars in both clusters is about 1 Myr. The age of evolved stars in Tr 16 is about 2.5 
Myr but we cannot estimate the upper limit of the age of Tr 14 because there is no evolved 
star in the cluster. 

Figure [11] shows the age distribution of PMS stars in the clusters. The average age of 
PMS members in both clusters are distributed between 1 Myr and 3 Myr. But, the age of 
the lower mass (logm < 0.2) stars seems to be older than 1 Myr, because the more luminous 
PMS stars (thus younger PMS stars) are preferentially detected due to our observational 
limit. On the other hand, the age of the massive PMS stars (logm ~ 0.3) is about 3 
Myr. These massiv e PMS stars in bo t h clu sters appear to still be in the Kelvin-Helmholtz 
contraction phase. 



Sung et al. 



(1997 



20041 ) suggested that the age of PMS star s in th e 



Hartmann 



(|2003|) 



Kelvin-Helmholtz contraction phase could be overestimated. Later 
discussed the age of PMS stars determined from PMS models. 

It is not easy to determine the age difference between these clusters because the age 
distribution of both clusters seems to be very similar. In order to estimate the difference in 
ages of the PMS stars in the clusters, deeper photometry down to sub-solar-mass stars is 
required. 



EDITOR: PLACE FIGURE [TT] HERE. 



4.5. Is Cr 232 an independent cluster? 

It is uncertain whether Collinder 232 (Cr 232), a poor young cluster with a small 
number of blue stars, located to the eastern side of Tr 14 is a separate cluster or an extended 
part of Tr 14. There are a few very early-type stars, including HD 93250 (03.5V((f+))) 

applied, HD 93162 is much closer to the 1 Myr isochrone in Figure 



-22 - 



and HD 303311 (05V), that are thought to be members of Cr 232 together with other 
proper motion members, but this simply suggests that the stars are at the same distance 
and have similar kinematic properties to Tr 14 and Tr 16. The existence of 03.5V((f+)) 
stars in Cr 232 (HD 93250) and Tr 14 (HD 93128) also implies that the two clusters 
are indistinguishable in age. Pre v ious s tudies have argued that Cr 232 may not be an 



independent cluster. iTapia et al. 



cluster (see Figure 13 of 



Tapia et al 



(2003) cone 



(|2003|)) 



uded that Cr 232 ap pears not to be a real 



Carraro et al 



( 120041 ) also suggested that Cr 



232 could be an extended part of Tr 14 rather than an independent cluster. We agree with 
these conclusions and have treated Cr 232 as an outer region of Tr 14 and Tr 16 rather than 
a separate cluster. 



5. Summary 

We have presented new UBVI C photometric data down to V = 19 for Tr 14 and Tr 16 
obtained with the lm telescope at Siding Spring Observatory and a SITe 2k CCD camera. 
We selected the members and c andidates using the X-ra y sour ce catalog from Chandra 



X-ray observations provided by lAlbacete-Co 
point-source catalog provided by 



Broos et al. 



ombo 



et al 



( 120081 ) and the CCCP X-ray 



(120111 ). the 2M ASS point-source catalo g, and 



the membership probability from the proper motion study of 



Cudworth et al 



(11993( 1 . 



We derived an abnormal reddening law in the rj Carina nebula and determined the 
foreground reddening to be E(B — V)f g = 0.36 ± 0.04 mag. An abnormal total-to-selective 
extinction ratio in the nebula, Rv,d = 4.4 ± 0.2, was determined using 141 early-type 
members. We derived the same distance modulus of Vq — My = 12.3 ±0.2 mag (d = 2.9 ±0.3 
kpc) for both Tr 14 and Tr 16, and concluded that the clusters are at the same distance. 



We derived the IMF and calculated the slope T = — 1.3 ± 0.1 for the rj Carina nebula 
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region down to our observational limit (logm ^ 0.2). The slope of the IMF, T, of the 77 
Carina nebula is very similar to that of NGC 6231 and is slightly shallower than that of 
NGC 2264, the Solar neighborhood, and the Pleiades. 

We estimated the age of MS stars in both clusters to be about 1 Myr from the 
comparison between the theoretical evolution models and the H-R diagram of the cluster. 
The more evolved star HD 93162 in Tr 16 is matched well with the isochrone of age 2.5 
Myr. We have also estimated the age of the PMS stars in the rj Carina nebula to be 1 - 3 
Myr. 

The authers thank the refree, D. Gies, for many insightful comments. This work is 
supported by a National Research Foundation of Korea (NRF) grant funded by the Korea 
Government (MEST) (Grant No. 20110114136). 
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Fig. 1. — Comparison of our photometry with four previous sets of CCD photometry. Eac h 



pan el from upper to low er shows t he difference be tween us and (1) 



(2) IVazquez et al. 



f)1996h . (3) iTapia et all (120031 ). and (4) 



Massey fc Johnson! ( 119931 ). 



Carraro et al. 



(120041 ). respectively. 
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Fig. 2. — 2MASS color-color diagra m of stars with ej J — H : H — Ks) < 0.1 ma g. The solid 



and d ashed lines are the MS relation (ISung et al. 



20081 ) and the reddening vector ( IFitzpatrick 



19991 ) for Ay = 10 mag. A triangle (red), filled circle (black), and cross (blue) represent NIR 



excess stars, proper motion members (P^70%) , and X-ray emission stars, respectively. 
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V-I B-V U-B 



Fig. 3. — Color-magnitude diagrams of stars with e ^ 0.1 mag. An open circle, filled circle, 
cross (blue), and triangle (red) represent data from photoelectric photometry, proper motion 
member, X-ray emission star, and NIR excess star, respectively. The stars with multiple 
membership criteria are marked with the symbol according to their primary membership 
criterion. The dashed and solid lines represent, respectively, the PMS locus and the ZAMS 
relation with Vq — My = 12.3 mag and E(B — V) = 0.61 mag. 
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B-V 

Fig. 4. — (U — B) vs. (B — V) color-color diagram of the stars with V ^ 17 mag. Open 
squares represent members selected from the reddening law (see Section H~Tj) . The solid line is 
the unreddened MS relation, while dotted lines are reddened MS relations with E(B — V) = 
0.36 and 0.9 mag. Other symbols are the same as Figure [31 The large open circle near 
E(B — V) = 0.9 mag represents rj Carinae. 
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Fig. 5. — Reddening law of early- type stars in the r/ Carina nebula. Solid lines correspond 
to Ry = 3.1 and Ry,d = 4.4 for E(B — V) ^ E(B — V)f g , and dashed line corresponds to 
E{B — V) = 0.9. Filled circles, open squares, and dots represent the stars with P M ^ 70%, 
the members selected from the reddening law (see Section I4TTT) . and the early-type stars with 
< 70% or without proper motion membership probability. 
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Fig. 6. — Smoothed reddening map of the observed region. The lines represent E(B — V) = 
0.5 (dashed line), 0.6 (thin solid line), and 0.7 (thick solid line), respectively. A R.A. and A 
decl. represent the angular distance in minutes of arc from 77 Carinae (R.A.= 10M5™3!591, 
decl= -59°41'4 / /26, J2000.0). 
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Fig. 7 .— Spatial distribution of member and candidate stars. Filled circles represent 
members from the proper motion study (P M ^ 70%), and from the abnormal reddening law 
in Figure |5j A R.A. and A decl. represent the angular distance in minutes of arc from rj 
Carinae. The members with X-ray emission or NIR excess are marked with a cross (blue) or 
a triangle (red). The stars with multiple membership criteria are marked by their primary 
membership criterion mentioned above. The size of symbols of members is proportional 
to the V magnitude of the stars. Small dots represent the candidate stars. The large 
circles represent the limit of Tr 14 (AR.A.= — 8.'1, Adecl.= 7 '.2, radius = 4'.2) and Tr 16 
(AR.A.= —1'.5, Adecl.= —2'.0, radius = 6f0). The rectangle to the west of Tr 14 represents 
the comparison field region used to remove the foreground contamination from the IMF (see 
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Fig. 8. 
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ZAMS fit t ing. O pen symbols repre sent t he known binaries f rom | Gamen et al. 



Mason et al 



(120091 ) 



Nelan et al 



Rauw et al. 



(12001 



20091 1. From the 



(12004] ) . and 

distribution of distance moduli of individual member stars, we obtain the same distance 
modulus (Vq — My = 12.3 mag) both for Tr 14 and Tr 16. The thick line represents the 
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Fig. 9. — H-R diagram of members of the two clusters. Dots, stars, filled circles, open circles, 
and open triangles represent PMS members, evolved stars (WN7, LBV), and MS members 
of Tr 14, Tr 16, and other area in th e observed region, res pectively. The thick and thin solid 



li nes represen t the Z AMS relation of 



Schaller et al. 



of 



Siess et al 



(119921 ) and the PMS evolutionary tracks 



(l2000h . The thick and thin dashed lines represent the isochrones interpolated 



from the PMS evolutionary tracks and the MS and post-MS evolutionary models. 
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Fig. 10. — IMF of Tr 14, Tr 16, and all stars in the observed region. Filled circles and open 
circles represent, respectively, the IMF for members and candidates (Table |3] case (2)), and 
for members-only (Table |3] case (1)). The error bars are based on >/~N. 
(i) The IMF calculated using the non-rotating models of 



Schaller et al. 



(ii) The IMF calculated using th e rotating models with the initial rotational velocity of 330 



km s 1 from 



Brott et al 



(j201lf ) for 5-6OM and the non-rotating models of 



Schaller et al. 



( 119921 ) for the other stars. 



(iii) The IMF calculated using the rotating models of 



Ekstrom et al. 



(1201 lh for all stars. 
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Fig. 11.— Age distribution of PMS members in Tr 14 and Tr 16. The size of dots is 
proportional to the number of PMS stars used in statistics. The dashed line represents the 
completeness limit of this work (logm = 0.2). For low-mass PMS stars, only young bright 
PMS members were observed and the age is therefore underestimated. For high-mass PMS 
stars, PMS evolutionary models overestimate the age (see Section |4~~4"1) . 



Table 1. Photometric Data 1 



ID a/2000 <5j2000 V U-B B-V V-I t v e v . B e fl _y ev-I N o6s P M X remark 2 2MASS 



1655 


10 


44 


00 


5 


-59 


33 


03.9 


14.775 


0.271 


0.464 


0.560 





007 


0.021 


0.012 


0.024 


6 


2 


6 


6 




X 


10440043-5933030 


1656 


10 


44 


00 


5 


-59 


32 


42.8 


18.538 






2.267 





026 






0.053 


4 








1 




X 


10440047-5932431 


1657 


10 


41 


00 


6 


-59 


48 


54.6 


16.546 




1.596 


2.435 





004 




0.009 


0.008 


6 





5 


6 






10440062-5948544 


1658 


10 


44 


00 


6 


-59 


34 


04.3 


13.834 


0.154 


0.483 


0.814 





002 


0.009 


0.006 


0.008 


6 


2 


6 


6 


19 


1 


10440061-5934043 


1659 


10 


41 


00 


6 


-59 


31 


52.3 


13.415 


-0.269 


0.287 


0.464 





007 


0.015 


0.011 


0.010 


6 


2 


6 


6 


58 




10440063-5931524 


1660 


10 


44 


00 


6 


-59 


32 


07.0 


17.758 




1.056 


1.147 





014 




0.024 


0.022 


4 





4 


4 








1661 


10 


41 


00 


7 


-59 


34 


17.3 


15.053 


0.467 


1.097 


1.496 





002 


0.037 


0.011 


0.004 


6 


1 


6 


6 


61 




10440067-5934173 


1662 


10 


44 


00 


7 


-59 


32 


33.5 


14.957 


0.177 


0.681 


0.740 





008 


0.046 


0.012 


0.020 


6 


2 


6 


6 






10440072-5932334 


1663 


10 


41 


00 


7 


-59 


49 


51.8 


12.495 


-0.283 


0.161 


0.257 





007 


0.016 


0.018 


0.008 


2 


2 


2 


2 






10440071-5949518 


1664 


10 


44 


00 


8 


-59 


45 


03.8 


13.951 


1.178 


1.320 


1.421 





002 


0.023 


0.005 


0.004 


6 


1 


6 


6 







10440079-5945038 


1665 


10 


44 


00 


8 


-59 


34 


38.9 


16.631 




1.519 


1.976 





014 




0.022 


0.019 


6 





5 


6 




X 


10440081-5934390 


1666 


10 


41 


00 


8 


-59 


37 


55.4 


18.387 




1.070 


1.726 





019 




0.071 


0.025 


4 





1 


4 






10440084-5937556 


1667 


10 


44 


00 


8 


-59 


48 


59.4 


15.895 


0.767 


0.661 


1.282 





Oil 


0.038 


0.013 


0.032 


6 


1 


6 


3 




D 


10440080-5948592 


1668 


10 


41 


00 


9 


-59 


33 


51.1 


17.605 




1.626 


2.194 





008 




0.058 


0.010 


4 





1 


4 




X 


10440091-5933512 


1669 


10 


11 


00 


9 


-59 


48 


01.1 


18.911 






1.657 





021 






0.026 


1 








1 






10440091-5948011 


1670 


10 


44 


00 


9 


-59 


35 


45.7 


10.670 


-0.625 


0.380 


0.712 





009 


0.012 


0.012 


0.012 


4 


1 


4 


4 


96 


X 


10440093-5935458 


1671 


10 


11 


00 


9 


-59 


17 


12.8 


16.048 


0.278 


0.825 


0.986 





004 


0.036 


0.007 


0.006 


6 


1 


6 


6 






10440093-5947126 



Note. 



Table[T]is published in its entirety in the electronic edition. A portion is shown here for guidance regarding its form and content. 



1 Membership probability (P M in percent) from proper motion study i Cudworth et al.lll993T ). 

2 D - The PSF of the object shows a convolution of two PSFs, but the object is measured as a single star for the sake of photometric error. X - X-ray emission 
stars. I - NIR excess stars. 



Table 2. Comparison of Photometry 



Paper 


AV 1 


N(m) 2 


A(B-V) 1 


N(m) 2 


A(CZ-B) 1 


N(m) 2 


A(V-Z) 1 


N(m) 2 


Feinstein ( 19631 


-0.006 ±0.124 


8(1) 


0.026 ±0.023 


8(1) 


0.034 ± 0.029 


8(1) 






Feinstein (1969) 


-0.006 ±0.128 


35(2) 


0.007 ±0.027 


32(5) 


0.014 ±0.034 


34(3) 






Feinstein et al. (1973) 


-0.041 ±0.072 


60(5) 


0.018 ±0.022 


58(7) 


-0.012 ± 0.041 


57(8) 






Herbst (1976) 


0.007 ± 0.064 


11(2) 


-0.006 ± 0.022 


12(1) 


-0.032 ± 0.046 


13(0) 






Feinstein (1982) 


-0.022 ±0.062 


26(5) 


0.028 ±0.035 


28(3) 


0.026 ± 0.061 


26(3) 


0.016 ±0.036 


38(12) 


Massev & Johnson ( 1993) 


-0.003 ±0.132 


160(23) 


0.031 ±0.114 


171(12) 


-0.033 ± 0.076 


127(56) 






Vazquez et al. (1996) 


-0.006 ± 0.022 


56(8) 


0.030 ±0.030 


55(9) 


0.018 ± 0.051 


57(6) 


0.007 ±0.031 


42(14) 


Tapia et al. (2003) 


-0.001 ±0.039 


345(44) 


-0.045 ± 0.072 


371(18) 


0.080 ± 0.087 


140(18) 


-0.042 ± 0.063 


356(33) 


Carraro et al. (2004) 


-0.018 ±0.030 


138(11) 


0.014 ±0.037 


140(9) 


0.060 ± 0.069 


134(14) 


0.072 ± 0.044 


177(18) 



O 

This - Others 



2 N and m represent number of compared stars and excluded stars. 
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Table 3. Comparison of the IMF slope 



Case Mass Range Tr 14 Tr 16 All Region Note 

(1) (logm^0.2) -1.3 ±0.1 -1.3 ±0.1 -1.4 ±0.1 Member + candidate 

(2) (logra^0.2) -1.3 ±0.1 -1.1 ±0.1 -1.2 ±0.1 Member only 

(3) (logm ^ 0.2) -1.3 ± 0.1 -1.2 ± 0.1 -1.3 ± 0.1 Member + candidate + field correction 



